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EXECUTIVE  SUMMARY 


Current  standoff  detection  protocols  for  the  nerve  agent,  O-ethyl-S-2- 
(diisopropylamino)ethyl  methylphosphonothiolate  (VX),  assume  that  the  material  is  of  high 
purity  and  that  the  speetra  are  dominated  by  VX  in  the  vapor  state.  We  report  the  acquisition  of 
spectra  of  weapons-grade  VX  in  the  mid-infrared  region  at  23  °C.  The  chemical  used  to  generate 
the  speetra  was  obtained  from  a  ton  container  that  was  analyzed  by  gas  chromatography-mass 
spectrometry  and  nuclear  magnetic  resonance  spectroscopy  to  determine  its  composition.  In 
contrast  to  the  pure  component  vapor-phase  reference  speetrum,  obtained  at  69  °C,  analysis  of 
the  data  from  the  weapons-grade  material  showed  that  the  speetra  were  dominated  by  the  volatile 
impurities  in  the  VX — primarily  acetone,  1.3-diisopropylcarbodiimide,  and  diisopropylanune. 
Concentrations  of  the  compounds  observ  ed  in  the  vapor  were  similar  to  the  predictions  made 
using  Raoult’s  law.  The  availability  of  the  absorptivity  coefficients  of  most  of  the  compounds 
that  were  found  in  the  effluent  enabled  us  to  calculate  a  mass  balance  that  encompassed  Sl%  of 
the  evaporated  vapor  from  the  saturator  cell  during  the  time  of  the  study. 
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PREFACE 


I  he  work  described  in  this  report  was  performed  under  the  direction  of  the 
Detection  Capability  Officer,  Defense  Threat  Reduction  Agency  Joint  Science  and  Technology 
Office  (DTRA  JSTO).  This  work  was  started  in  February  2009  and  completed  in  March  2009. 

The  use  of  either  trade  or  manufacturers’  names  in  this  report  does  not  constitute 
an  official  endorsement  of  any  commercial  products.  This  report  may  not  be  cited  for  purposes 
of  advertisement. 

This  report  has  been  approved  for  public  release.  Registered  users  should  request 
additional  copies  from  the  Defense  Technical  Information  Center;  unregistered  users  should 
direct  such  requests  to  the  National  Technical  Information  Service. 
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VAPOR-PHASE  INFRARED  SPECTRAL  STUDY  OF  WEAPONS-GRADE 
O-ETHYL  S-2-(DIISOPROPYLAMlNO)ETHYL 
METHYLPHOSPHONOTHIOLATE  (VX) 


1.  INTRODUCTION 

Among  the  ehemieal  warfare  agents,  the  organophosphonates  (classified  as  nerve 
agents)  eneompass  a  elass  of  highly  toxie  compounds  w  ith  a  wide  range  of  persistencies  and 
toxieities.  The  nerve  agents  ean  be  further  split  into  three  subclasses  with  similar  structures: 

(1)  O-alkyl  V,V-dialkylphosphoramidoeyanidate  (e.g.,  GA  or  tabun);  (2)  alkyl 
alkylphosphonofluoridate  (e.g.,  GB  or  sarin);  or  (3)  0-alkyl-S-2-(dialkylaminoethyl) 
alkylphosphonothioate  (e.g.,  the  V-series). 

This  report  eovers  a  study  of  the  infrared  (IR)  speetra  produced  by  evaporation  of 
a  member  of  the  third  subelass  listed  in  the  paragraph  above,  ()-ethyl-S-(2- 
diisopropylamino)ethyl  methylphosphonothiolate  (VX).  VX  has  the  Chemieal  Abstracts  Serviee 
(CAS)  Registry  Number  57082-69-9.  The  molecular  formula  is  Ci  i  H26NO2PS,  and  the 
molecular  weight  is  267.37.  The  structure  is  shown  in  Figure  I . 


Figure  1.  Structure  of  ()-ethyl-S-(2-diisopropylamino)ethyl  methylphosphonothiolate  (VX). 

Standoff  detection  systems,  especially  those  operating  in  the  mid-IR  region  of  the 
spectra,  are  primarily  programmed  to  identify  ehemieal  warfare  agents  on  the  basis  of  their 
vapor-phase  1R  signatures.  Several  of  the  G-series  agents  (subclasses  1  and  2  above)  have 
significant  vapor  pressures  ( P )  and  volatilities  (GB:  Pisx-  2.48  torr1)  that  facilitate  the 
detection  of  a  vapor  eloud.  However,  the  vapor  pressures  and  volatilities  of  the  subelass  of  V 
agents  tend  to  be  lower  than  those  of  the  G  agents.  VX  has  a  vapor  pressure  of  0.00088  torr  at 
25  °C,2  whieh  is  more  than  2800  times  lower  than  sarin.  Furthermore,  the  vapor  pressures  of 
many  of  the  impurities  in  VX  are  mueh  higher  than  that  of  the  pure  agent.  A  first-order 
approximation  would  prediet  that  impurities  in  VX,  even  at  very  low  concentrations,  eould 
contribute  to  the  IR  signature  in  a  manner  that  is  disproportionate  to  their  molar  or  mass 
concentrations.  Therefore,  the  compound  elearly  presents  a  challenge  for  standoff  detection. 
Nevertheless,  with  a  lethal  dose  for  50%  of  the  population  (LD50)  of  5  mg,1  the  early  warning  of 
potential  contamination  is  eritieally  important. 
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2. 


EXPERIMENTAL  PROCEDURES 


2.1  Instruments  and  Laboratory  Apparatus 

The  system  used  to  generate  the  continuous  vapor  stream  was  an  adaptation  of  the 
saturator  cell  method  developed  at  the  U.S.  Army  Edgewood  Chemical  Biological  Center 
(ECBC)  for  measuring  the  volatility  of  chemical  warfare  (CW)  agent  related  compounds.1 

The  method,  which  was  modified  to  generate  continuous  streams  of  chemical 
compounds  and  obtain  quantitative  vapor-phase  1R  spectra,  has  been  used  to  measure  the 
absorptivity  coefficients  of  a  variety  of  chemical  agent  related  compounds.  The  saturator 
passes  a  stream  of  nitrogen  carrier  gas,  which  was  obtained  from  the  boiloff  of  a  bulk  liquid 
nitrogen  tank,  across  an  alumina  Soxh let-shaped  wick  positioned  in  a  glass  holder  filled  with  the 
analyte.  This  technique  yields  a  saturated  vapor-liquid  equilibrium  of  the  analyte  on  the 
downstream  side  of  the  saturator  cell  at  flow  rates  that  arc  less  than  approximately  100  mL/min. 
When  the  saturator  cell  is  suspended  in  a  constant-temperature  bath,  the  concentration  of  the 
analyte  can  be  predicted  by  its  vapor  pressure  at  the  temperature  of  the  bath.  At  carrier  rates 
>100  mL/min,  the  efficiency  of  the  saturator  cell  will  slowly  decline. 

Using  the  saturator  cell  and  bath,  we  have  been  able  to  calculate  the  change  in  the 
vapor-liquid  equilibrium  (percent  of  saturation)  as  a  function  of  the  carrier  rate  (Figure  2)  by 
measuring  the  change  in  mass  at  a  range  of  carrier  rates.  The  mass  rate  will  remain  constant  for 
a  pure  compound  (assuming  that  the  carrier  rate  and  temperature  do  not  change).  The  carrier  gas 
was  controlled  with  a  Brooks  Model  5850S  (Emerson,  Hatfield.  PA)  mass  flow  controller.  A 
second  mass  flow  controller  was  used  to  add  diluent  to  the  stream,  which  provided  an  additional 
way  to  adjust  the  concentration  of  the  compound  delivered  to  the  White  cell  of  the  Fourier 
transform  infrared  (FTIR)  spectrometer.  The  linearity  of  the  S-series  mass  flow  controllers  was 
adjusted  using  a  second-order  polynomial,  which  resulted  in  rate  uncertainties  of  approximately 
1%  or  better  at  flows  >25%  of  full  scale. 

Spectra  were  obtained  with  a  Brukcr  model  1FS/66V  FTIR.  The  instrument 
was  equipped  with  deuterated  triglyeinc  sulfide  and  mercury  cadmium-telluride  (HgCdTc) 
detectors  and  was  capable  of  obtaining  spectra  with  a  maximum  spectral  resolution  of 
0.1 125  cm  1  (unapodized).  The  interferograms  were  recorded  from  (15,798  to  0)  cm  1  with  a 
resolution  of  0.5  cm  l.  Single-beam  spectra  were  averaged  from  128  scans.  A  background 
spectrum  of  the  dry  nitrogen  was  taken,  the  carrier  was  started,  the  vapor  from  the  saturator  cell 
was  switched  to  the  White  cell,  and  the  acquisition  of  sample  spectra  was  started  at  intervals  of 
3  min  until  the  carrier  flow  was  stopped.  Absorbance  (logio)  spectra  were  processed  with  boxcar 
apodization  and  2x  zero  filled  to  obtain  a  data  spacing  of  0.25  cm  The  instrument  is  equipped 

with  a  variable  path  White  cell.  The  experimental  data  used  a  path  length  of  2.727  m.  The 
temperature  of  the  White  cell  was  maintained  at  23  ±  0.1  °C  through  the  use  of  a  thermostatically 
controlled  chamber  that  enclosed  the  spectrometer  and  cell.  Data  were  acquired  at  a  speed  of 
60  kHz  (helium-neon  [HeNe]  laser  zero-crossing  frequency)  using  the  HgCdTe  detector. 
Absorbance  spectra  of  the  vapor  effluent  were  computed  using  background  spectra  of  clean,  dry 
nitrogen.  To  minimize  the  effects  of  nonlinearity  in  the  detector,  the  interferograms  were 
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processed  using  the  proprietary  Opus  nonlinearity  correction  function.  All  interferograms  were 
archived,  thus  enabling  additional  post-processing  of  data  as  needed. 


Figure  2.  Saturator  cell  efficiency  as  a  function  of  carrier  rate.  The  equation  for  the  fitted  line 
above  is  E  =  (0.0001 402  *  R)  +  1 .012,  where  E  =  fractional  efficiency  and  R  =  carrier  rate 
referenced  to  a  temperature  of  23  °C. 


The  saturator  cell  was  weighed  prior  to  the  start  of  each  experiment.  At  the  end 
of  each  day’s  experiment,  the  saturator  cell  was  dried,  equilibrated  for  a  minimum  of  2  h  at  room 
temperature,  and  then  reweighed.  The  change  in  mass  was  compared  to  the  data  obtained  from 
the  IR  spectra  of  the  effluent  (described  in  Section  2.3). 

Temperature  and  pressure  data  were  recorded  using  National  Institute  of 
Standards  and  Technology  (NIST)-traceable  digital  manometers  and  thermometers,  and  all  data 
were  archived. 

2.2  Analysis  of  Feedstock 


The  material  used  to  generate  the  vapor  streams  for  the  experiments  was  an 
archived  sample  from  35274-15,  Ton  Container  A49.  An  aliquot  of  the  sample  was  analyzed  by 
gas  chromatography-mass  spectrometry  (GC-MS)  in  split  mode  on  13  September  2007. ^  Results 
are  shown  in  Table  1 .  Fractions  shown  in  the  table  were  computed  from  the  areas  in  the  GC 
chromatogram.  Although  these  can  provide  an  indication  of  the  relative  concentrations  of  the 
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compounds,  the  fractions  are  only  an  approximation  of  either  mass  or  molar  fractions  of  the 
materials. 


Analysis  of  the  vapor-phase  1R  spectra  of  the  effluent  from  the  VX  showed  the 
presence  of  an  additional  compound,  acetone  (which  was  not  included  in  Table  1).  The  original 
source  of  the  acetone  in  the  sample  was  not  determined,  although  suggested  sources  included 
residue  from  the  use  of  acetone  to  rinse  the  ton  containers  or  the  glassware  that  was  used  to  take 
the  samples  from  the  ton  containers.  Analysis  by  nuclear  magnetic  spectroscopy  (NMR.  1 H) 
confirmed  the  identification  at  a  mass  fraction  of  0.27  ±  0. 12%  (mol  fraction  1.23  ±  0.55%,  with 
a  confidence  interval  of  approximately  95%).  When  the  results  from  the  GC-MS  analysis  were 
reexamined  using  an  extracted  ion,  the  peak  for  acetone  was  found;  although  the  peak  was  too 
small  to  determine  the  area  fraction. 


Table  1.  Summary  of  Results  from  GC-MS  Analysis  of  VX  Sample  from  Ton  Container  A49. 
Examination  of  the  peak  area  fractions  from  the  chromatogram  gives  only  an  approximation  of 
mass  or  molar  fractions  of  the  individual  compounds  found. 


Reaction 
Time  (RT) 
(min) 

Molecular 

Weight 

Compound 

Peak  Area 
(%) 

2.38 

101 

Diisopropylamine  (D1A) 

0.31 

7.83 

126 

1,3-Diisopropylcarbodiimide  (D1CD1) 

0.84 

11.13 

168 

Diethyl  methylphosphonate  (DEMP) 

0.10 

12.17 

161 

2-(Diisopropylamino)cthanethiol 

0.32 

13.42 

187 

2-(Diisopropylamino)ethyl  vinyl  sulfide 
[(iPr)2-N-CH2CH2SCH=CH2l 

0.07 

13.65 

144 

AW’-Diisopropylurea 

0.51 

14.73 

189 

2-Diisopropylaminoethyl  ethyl  sulfide 

0.05 

16.29 

230 

Diethyl  dimethypyrophosphonate 

1.14 

18.40 

251 

O-Ethyl  0-[2-(Diisopropylamino)ethyl] 
methyl  phosphoric  acid 

0.05 

20.09 

267 

VX 

93.55 

20.37 

- 

Unknown  (base  1 14  m/z) 

1.87 

21.46 

281 

O-Ethyl  S-[2-(diisopropylamino)ethyl] 
cthylphosphonthiolate 

0.40 

23.65 

320 

Bis[(diisopropylamino)ethyl]  disulfide 

0.79 

14 


2.3 


Theory  and  Calculations 


The  pressure  over  a  solution  is  equal  to  the  total  partial  pressures  of  the  solutes  in 
the  solution.  According  to  Raoult’s  law,  the  partial  pressure  of  a  solute  is  proportional  to  the 
partial  pressure  of  the  pure  solute  and  its  molar  fraction  in  the  solution:1" 


P  =  PaXa+PbXb+P(Xc...  (I) 

where  P  is  the  pressure  over  the  solution,  and  PaXa...  arc  the  pressure  of  the  pure  solvent  and  its 
mole  fraction.  For  dilute  solutions  with  only  two  components,  particularly  when  the  two 
compounds  are  chemically  similar,  the  use  of  Raoult’s  law  can  yield  predicted  pressures  that  arc 
close  to  the  actual  values.  In  complex  solutions,  with  a  variety  of  intermolccular  interactions,  or 
for  azeotropic  solutions,  pressure  can  only  be  approximated  by  applying  Raoult’s  law  and,  in 
general,  it  needs  to  be  determined  empirically.  The  equation  can,  however,  be  used  as  a  starting 
point  for  estimating  the  evaporation  of  a  solution  and,  in  this  case,  the  solution  was  weapons- 
grade  VX  from  a  saturator  cell. 

As  stated  previously,  saturator  cells  have  been  used  at  ECBC  to  measure  the 
vapor  pressure  of  a  variety  of  compounds,  including  CW  agents  in  the  solid  phase.1 1  Eq  2  below 
is  used  to  calculate  the  vapor  pressure  of  a  compound  ( Pa )  from  saturator  cell  measurements 


P=P. 


n. 


"earner 


(2) 


where  P is  the  total  pressure  in  the  saturator  cell,  na  is  the  number  of  moles  of  the  compound 
(from  the  change  in  mass  and  molecular  weight  of  the  compound),  and  wcarTjCT  is  the  number 
moles  of  carrier  gas  (obtained  from  the  flow  rate  of  the  carrier  and  the  ideal  gas  law) 


FRAt-P 


RT 


(3) 


where  FR  is  the  flow  rate  of  the  carrier  in  L  min  ',  t=  time  in  minutes,  P  =  101 ,325  Pa, 

R  —  8,3 14.5  L  •  Pa  (mol  •  K)  *,  and  T  =  the  calibration  temperature  of  the  mass  flow  controller  in 
Kelvin  (generally  294.26  K  for  the  mass  flow  controllers  in  our  laboratory). 


For  low-volatility  compounds,  where  fhamcJna  approaches  1,  na  in  the  lower  part 
of  the  fraction  in  eq  2  can  be  ignored,  and  the  equation  can  then  be  rearranged  and  combined 
with  eq  1  to  yield 


P.Xn  . 


_  a  a  earner 

°  P 

1  sal 


(4) 


By  setting  t  =  I  in  eq  3,  as  long  as  Pa  is  known,  and  Xa  =  l  (i.e.,  the  purity  of  the  compound  in  the 
saturator  cell  is  100%),  eq  4  can  be  used  to  predict  the  evaporation  rate  of  a  compound  from  a 
saturator  cell.  This  should  be  a  steady  state  as  long  as  PM  remains  constant.  In  most  cases, 
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saturator  cells  are  used  to  generate  vapor  from  compounds  having  purities  of  >95%.  In  such 
cases,  the  concentrations  of  other  compounds  in  the  solution  (A^)  are  typically  1  to  2%,  and  eq  4 
is  a  Raoult’s  law  prediction  of  the  instantaneous  evaporation  rate  from  a  saturator  cell.  In  most 
cases,  the  vapor  pressures  of  the  compounds  in  the  solution  will  differ.  At  least  for  an  ideal 
Raoult  solution,  examining  eq  4  suggests  that,  regardless  of  the  starting  molar  concentrations,  if 
PaXa  <  PfrXf,,  Xb  will  decrease  and  Xa  will  increase  over  time.  As  Xb  decreases,  its  evaporation 
rate,  «/,,  will  also  decrease  asymptotically.  In  the  simplest  case,  /(/)  =  ln(«f,),  whereof/)  can  be 
described  using  the  relationship  at  +  b,  which  can  then  be  written  as 

b0t +  bi  =\n(nh)  (5) 

Eq  5  can  then  be  rewritten  as 

cv+Al  =  nh  (6) 


The  integral  of  eq  6  is  computed  as 


,V+*i 


(7) 


The  evaporation  rate  of  a  compound  («/,)  can  be  determined  from  its  concentration 
in  the  vapor,  which  can  be  calculated  from  the  IR  spectra  by  rearranging  Beer’s  law12 


A,(F)  =  al(F)Lci 


(8) 


to  obtain 


4(P) 

a,(v)L 


(9) 


where  c,-  is  the  concentration  of  compound  /  in  the  vapor,  /f  (r)  is  the  absorbance  of  the 
compound  at  a  wave  number,  a (v)  is  the  absorptivity  coefficient  of  the  compound  at  that  wave 
number,  and  L  is  the  path  length. 


Eq  9  can  be  computed  using  A,(v)  and  a*  as  either  intensity  at  a  single  wave 

number  (peak  height)  or  integrated  across  a  region.  The  reference  spectra  were  recorded  at  a 
higher  resolution  (0. 125  cm  )  than  the  sample  spectra.  In  such  a  case,  integrating  across  a 
spectral  region  can  improve  the  accuracy  of  the  calculations.  For  that  reason,  we  preferred  to  use 
integrated  absorbance  to  calculate  the  concentrations  of  the  compounds.  The  exceptions  were 
(1)  Acetone:  An  unknown  compound  was  present  in  the  VX  with  an  absorption  band  adjacent  to 
the  ketone  carbonyl  stretch  band  of  acetone.  The  unknown  compound  apparently  had  a  low 
vapor  pressure,  and  its  intensity  changed  little  over  the  course  of  the  experiments.  Its 


*  This  also  implies  that  the  concentration  ofXa  and  proportion  of  rta  increase  gradually.  Because  we  arc  considering 
the  more  typical  starting  case  in  which  Xa  is  much  larger  than  Xh,  the  proportional  effect  on  the  rate  of  na  is  much 
smaller. 
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concentration  was  apparently  low,  and  it  did  not  interfere  with  the  acetone  except  at  lower 
concentrations  of  aeetonc,  which  were  then  calculated  using  peak  height.  (2)  DEMP:  The 
phosphoryl  band  used  to  calculate  the  concentration  of  the  compound  is  in  a  busy  region  of  the 
spectrum,  with  bands  from  several  other  compounds  nearby.  Concentrations  of  the  compound 
were  computed  using  peak  heights  only.  When  peak  heights  were  used,  they  were  computed  as 
means  across  a  narrow  region  (~1  to  2  cm  '),  to  minimize  the  effect  of  noise. 

The  reference  spectra  used  to  calculate  the  concentrations  of  compounds  in  the 
vapor  are  provided  in  units  of  reciprocal  micromoles  per  mole  x  meter  (pmol/mol)  'm  1  (parts 
per  million-meter,  ppm-m),  where  micromoles  =  MW  x  10  6  (yielding  a  mass  in  micrograms) 
and  moles  =  molar  volume  of  the  diluent  gas  (nitrogen  =  24.13  L  at  294.26  K,  which  is  the 
calibration  temperature  of  the  mass  flow  controllers).  The  reference  spectra  were  normalized  to 
296.15  K  (the  temperature  of  the  White  cell  used  in  our  measurements)  and  101,325  Pa.  The 
following  equation  can  thus  be  used  to  calculate  the  mass  rate,  R ,  of  <?,: 

R  c,{FR)(MW){P^) 

24.13(101325)  W 

where  R  =  instantaneous  mass  rate  in  micrograms  per  minute,  FR  =  the  total  flow  rate  of  carrier 
and  diluent  in  liters  per  minute  at  294.26  K  and  101,325  Pa,  MW=  molecular  weight  of  the 
compound,  and  Pamb  =  ambient  pressure  in  pascals.  Because  acquiring  128  scans  at  0.5  cm 
resolution  required  approximately  1  min,  R  was  actually  a  near-instantaneous  rate.  For  those 
compounds  having  quantitative  reference  spectra,  the  least-squares  fit  of  cq  5  was  then 
calculated,  substituting  R  for  n.  This  yielded  the  coefficients  bo  and  bt  for  /(/)  (eqs  6  and  7).  If 
J[t)  could  be  determined  for  each  compound  in  the  vapor  from  the  saturator  cell,  then  this  should 

permit  the  calculation  of:  II  m,  ,  where  m,  is  the  total  mass  of  each  compound  evaporated  from 

the  saturator  cell.  In  other  words,  we  believed  that  it  might  be  possible  to  compute  a  mass 
balance  from  the  1R  spectra  that  could  then  be  compared  to  the  gravimetric  data  that  was 
obtained  by  weighing  the  saturator  cell. 

The  data  ultimately  encompassed  500  spectra,  from  12  experiments  that  were 
conducted  over  a  total  of  more  than  27  h,  from  which  5  compounds  were  identified  and 
computed.  This  clearly  could  have  been  an  overwhelming  computational  task.  A  customized 
algorithm  was  written  in  MatLab  (Mathworks,  Natick,  MA)  to  perform  batch  processing  of  a 
number  of  tasks  (baseline  corrections  of  spectra,  spectral  integration,  and  Beer’s  law 
computations  of  concentrations),  thus  greatly  simplifying  the  computations. 


3.  RESULTS  AND  DISCUSSION 

3. 1  Prediction  of  Initial  Composite  Spectrum 

Prior  to  running  experiments  with  actual  VX,  we  attempted  to  predict  a  composite 
spectrum  using  Raoult’s  law  and  referring  to  a  report  that  gave  the  mean  concentrations  of 
compounds  found  in  60  ton  containers  of  the  agent.1  Those  reported  concentrations  were 
generally  similar  to  those  in  Tabic  1.  Among  the  volatile  compounds  having  available  reference 
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spectra,  the  concentrations  of  two  compounds  were  lower:  D1A  (0. 14%)  and  DEMP  (0.06%): 
and  the  concentrations  of  two  other  compounds  were  higher:  DICD1  (1 .74%)  and 
2-(diisopropylamino)ethanethiol  (0.89%).  Acetone  was  reported  at  a  mean  concentration  of 


0.01%. 


The  absorptivity  coefficient  (a)  of  VX  was  obtained  from  the  database  compiled 
by  Pacific  Northwest  National  Laboratory  (PNNL).14  Absorptivity  coefficients  of  other 
compounds  were  obtained  from  data  that  were  acquired  in  our  laboratory.1  Vapor  pressures 
were  computed  from  a  NIST  vapor  pressure  database  for  acetone  and  D1A1"  and  from  ECBC 
data  for  DEMP  and  VX/’  The  vapor  pressure  of  2-(diisopropylamino)ethanethiol  was 
computed  using  Antoine  coefficients  calculated  from  multiple  sources  of  data1K  22  as  well  as 
vapor  pressure  data  acquired  while  measuring  the  absorptivity  coefficient  of  the  compound  in 
our  laboratory.  The  vapor  pressure  of  D1CD1  was  calculated  from  coefficients  that  were 
computed  from  measurements  of  the  compound.  These  measurements  were  made  by:  ( 1 ) 
Chemical  Sciences  Division  at  ECBC‘  and  (2)  our  laboratory.  The  ECBC  Chemical  Sciences 
division  provided  vapor  pressure  measurements  for  D1CD1  that  were  similar  to  values  reported 
by  American  Cyanamid  Company  (in  the  temperature  range  of  51  to  146  °C  in  1962).24  Our 
measurement  of  the  vapor  pressure  of  the  compound  at  1 5  °C  was  concurrent  with  the 
acquisition  of  the  absorptivity  coefficient.1^  Table  2  provides  the  vapor  pressures  of  pure 
compounds,  at  the  saturator  cell  temperatures  from  this  work,  along  with  the  associated  data 
sources.  Figure  3  provides  the  composite  spectrum,  generated  from  the  Raoult’s  law  predictions 
of  the  vapor  pressures  of  the  volatile  compounds  (as  well  as  VX).  Figure  4  supplies  the  vapor- 
phase  spectrum  of  VX  from  the  PNNL  database,  acquired  at  69  °C.  Although  the  impurities 
were  present  at  low  concentrations,  Raoult’s  law  predicted  that  their  1R  signatures  would  be 
observable  in  the  vapor. 


Table  2.  Calculated  Vapor  Pressures  at  Temperatures  of  Selected  Compounds  Observed  in 
Vapor  from  Wcapons-Grade  VX.  The  pressures  provided  are  for  the  pure  compounds  at  the 


Compound 

Vapor  Pressure/Pa 

Reference 

5  °C 

15  °C 

23  °C 

Acetone 

12,088 

16 

D1A 

3906 

6613 

9772 

16 

D1CD1 

153 

317 

543 

15,  24 

DEMP 

10.8 

25.4 

48.0 

17 

2-(Diisopropylamino)ethanethiol 

12.1 

28.3 

53.5 

19-23 

VX 

0.011 

0.038 

0.094 

2 
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Figure  3.  Composite  spectrum  of  VX  generated  by  applying  Raoult’s  law  predictions  of  the 
vapor  pressures  and  absorptivity  coefficients  of  compounds  found  in  weapons-grade  VX. 
Concentrations  of  compounds  were  mean  values  for  ton  containers  of  VX  as  reported  in 
literature  reference  13. 


Figure  4.  Vapor-phase  absorptivity  coefficient  of  VX  from  PNNL  database,  valid  at  69  °C. 

Little  or  no  evidence  of  VX  is  shown  in  the  predicted  spectrum  of  the  vapor  from  weapons-grade 
VX  (Figure  3). 
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3.2 


Summary  of  Experimental  Conditions 


A  total  of  12  experiments  were  conducted  between  25  February  and  19  March 
2009.  Bath  temperatures  and  carrier/diluent  flow  rates  are  summarized  in  Table  3. 


Table  3.  Summary  of  Bath  Temperatures  and  Carrier/Diluent  Temperatures  in  the  VX  Vapor 
Experiments. 


Date(s) 

Bath  Temperature 
(°C) 

Flow  Rate 

(cm*  min  1  [seem  at  21.1  °C]) 

Carrier 

Diluent 

25  Feb  2009 

5 

50 

2500 

26  Feb  2009 

15 

50 

2500 

4-6,  10-12,  16  Mar 

23 

100 

2500 

17  Mar 

23 

500 

2500 

18  Mar 

23 

1000 

2500 

19  Mar 

23 

2000 

2000 

3.3  Results  from  the  Initial  Experiment 

The  first  day’s  experiment  was  run  at  a  bath  temperature  of  5  °C,  at  a  carrier  rate 
of  50  seem  (standard  cubic  centimeters  per  minute,  referenced  to  294.26  K.),  and  with  a  diluent 
of  2500  seem  to  capture  the  spectra  of  the  higher  volatility  compounds  present  in  VX,  which 
otherwise  might  have  purged  quickly  from  the  saturator  cell.  A  spectrum  of  the  effluent 
(acquired  at  15  min  into  the  experiment),  along  with  the  reference  spectra  of  acetone  and  DICDI 
for  comparison,  are  provided  in  Figure  5.  The  prominent  absorption  feature  associated  with  the 
ketone  carbonyl  (C=0  stretch)  from  acetone  was  unexpected,  given  the  initial  report  from  the 
GC-MS  analysis  of  the  VX  sample/  Furthermore,  evaluation  of  the  intensity  of  the  band 
indicated  that  the  acetone  was  likely  present  at  a  concentration  greater  than  the  mean  value  of 
0.01%  reported  in  literature  reference  13,  which  contained  data  that  was  likely  to  have  been 
obtained  by  GC-MS  analysis  also. 

A  bias  in  the  concentration  of  the  acetone  reported  in  the  literature  could  have 
occurred  for  two  reasons:  (I)  The  numbers  that  were  obtained  during  integration  of  a 
chromatogram  from  GC-MS  analysis,  represented  area  fractions  of  the  peaks  from  the  individual 
compounds  in  the  sample.  The  peak  area  fractions,  in  turn,  were  calculated  from  the  totals  of  the 
charged  mass  fragments  detected.  They  were  not  a  direct  measurement  of  the  portion  of  either 
the  mass  (or  molar  concentration)  of  the  compounds  making  up  the  sample.  Compounds  that 
fragment  similarly  to  one  another  in  the  electron  beam  may  have  peak  area  fractions  that  are 
similar  to  their  respective  molar  concentrations,  although  this  is  not  guaranteed.  (2)  To  detect 
the  presence  of  the  compounds  with  the  lowest  boiling  temperatures,  GC  analysis  must  be 
performed  with  the  neat  sample  undiluted  in  solvent.  This  requires  the  GC  to  be  operated  in 
“split”  mode,  in  which  most  of  the  sample  that  is  injected  into  the  heated  liner  in  front  of  the  GC 
column  is  vented  to  bypass  the  column.  The  low  boiling  temperature  of  acetone  (56.0  °C)  may 
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have  eaused  mueh  of  the  chemical  to  vaporize  and  expand  into  the  liner  more  quiekly  than  the 
other  compounds  in  the  VX,  causing  most  of  it  to  be  vented. 


x  io'5 


Figure  5.  Spectrum  of  effluent  from  the  VX-filled  saturator  eell  15  min  after  the  start  of  the 
experiment.  The  speetrum  is  dominated  by  DICDI  (2130  em  *)  and  aeetone  (1740  em  ').  The 
reference  spectra  of  the  two  compounds  have  been  offset  for  display  purposes. 

Acetone  at  a  mass  or  molar  fraction  of  0.01%  would  not  have  been  measurable  by 
NMR.  Evaluation  of  the  intensity  of  the  absorption  features  shown  in  Figure  5  indicated  that  it 
might,  nevertheless,  have  been  possible  to  detect  and  measure  acetone  with  NMR.  This  was 
subsequently  done  using  an  internal  standard  method.  The  reported  molar  fraction  of  0.27  ± 

0. 12%  for  acetone  has  a  large  uncertainty  because  of  the  relatively  small  size  of  the  acetone  peak 
in  the  NMR  spectrum  of  the  VX.  Nevertheless,  this  value,  combined  with  the  results  from  GC- 
MS  analysis,9  provided  a  new  starting  point  to  predict  the  speetrum  of  the  VX. 

Upon  analysis,  the  result,  shown  in  Figure  6,  still  appeared  to  overestimate  the 
relative  contribution  of  the  aeetone.  There  were  several  possible  explanations  for  this:  (1)  The 
actual  concentration  of  aeetone  eould  have  been  on  the  low  side  of  the  uncertainty  range  reported 
after  using  the  NMR.  (2)  The  vapor  pressure  of  acetone  in  solution  with  VX  may  not  have 
followed  an  ideal  Raoult’s  law  behavior.  (3)  The  molar  concentration  of  DICDI  was  higher  than 
the  peak  area  fraction  reported  as  a  result  of  using  the  GC-MS.  Unfortunately,  beeause  the 
protons  in  DICDI  are  similar  to  those  found  in  VX  and  other  nitrogen-eontaining  impurities,  it 
was  not  possible  to  obtain  an  independent  estimate  of  the  VX  concentration  using  NMR. 
Therefore,  the  true  mass  or  molar  concentration  of  the  stabilizer  in  the  VX  remained 
undetermined.  The  experience  did,  however,  illustrate  the  usefulness  of  a  multidisciplinary 
approach,  in  which  several  analytical  techniques  were  used  to  provide  pieees  of  the  puzzle  and 
emphasize  the  need  for  empirical  verifications  of  the  predictions. 
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Figure  6.  Composite  spectrum  of  VX  generated  by  applying  Raoult’s  law  predictions  of  the 
vapor  pressures  and  absorptivity  coefficients  of  compounds  found  in  weapons-grade  VX.  The 
spectrum  reflects  the  concentration  of  acetone  observed  by  the  use  of  the  NMR  and  the 
concentrations  of  other  compounds  observed  during  GC-MS  analysis  (literature  reference  9). 

3.4  Quantitation  of  Compounds  by  Vapor-Phase  IR:  Results  from 

the  Experiments  at  5  and  1 5  °C 

When  plotted  as  a  function  of  time,  the  mass  rate  of  acetone  obtained  during  the 
first  day’s  experiment  appeared  to  exhibit  the  decay  that  was  predicted  in  eq  5,  as  shown  in 
Figure  7.  The  calculated  mass  rate  of  acetone  at  to  =  50.7  pg-min  !  decayed  to  24.2  pg-min  1  at  / 
=  69.5  min,  which  was  the  end  of  the  experiment.  The  total  mass  of  acetone  observed  (from 
eq  7,  the  integral  of  eq  6)  was  2.5  mg.  The  mass  rate  of  D1CD1  remained  nearly  constant  at  a 
mean  rate  of  8.8  pg-min  ',  resulting  in  a  total  mass  of  0.6  mg.  Given  the  peak  area  fraction  of 
0.83%  for  the  compound  and  the  small  R  value,  it  would  have  been  difficult  to  detect  its  decay. 
The  total  decay  for  both  D1CDI  and  acetone  was  3.1  mg.  The  change  in  mass  obtained 
gravimetrieally,  by  weighing  the  saturator  cell  at  the  beginning  and  end  of  the  experiment,  was 
also  3.1  ±  0.2  mg.  The  initial  mass  rate  of  DIA  was  predicted  to  be  25  pg-min  1  when  Raoult’s 
law  was  applied.  (This  prediction  was  made  on  the  assumption  that  the  GC-MS  area  fraction  of 
0.31%  equaled  its  molar  concentration.)  The  absorption  feature  of  DIA  (centered  near  690  cm  ') 
that  was  used  for  quantitation  was  weak  [a  =  0.00021  (pmol/mol)  'm  ',  which  gave  a  prediction 
of  A  =  0.001  ]  and  close  to  the  detector  cutoff  (root  mean  square  [RMS]  noise  =  0.0002  A).  This 
was  consistent  with  the  failure  to  observe  DIA  in  the  vapor.  Because  the  prediction  obtained 
using  Raoult’s  law  for  the  mass  rate  of  VX  was  R  =  0.056  pg-min  ',  and  a  similar  prediction  was 
obtained  for  DEMP  (R  =  0.061  pg-min  '),  the  failure  to  detect  either  of  these  two  compounds 
was  also  expected. 
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Figure  7.  Mass  rate  of  aeetone  observed  in  effluent  from  VX-filled  saturator  eell  on 
25  February  2009.  The  deeay  in  aeetone  output  appeared  to  follow  the  behavior  predieted  using 
eq  5:  ln(7?)  =  hot  +  b\,  where  b0  =  -0.0165  and  b\  =  3.925. 

During  the  seeond  day’s  experiment,  the  constant  temperature  bath  was  set  to 
15  °C,  the  earrier  flow  was  operated  at  50  seem,  and  the  diluent  was  set  to  2500  seem.  The 
relative  intensities  of  aeetone  and  D1CDI  were  similar  to  those  observed  at  5  °C.  The  integrated 
masses  of  aeetone  and  DICD1  were  2.4  and  1.8  mg,  respectively.  In  contrast  to  the  previous 
day’s  experiment,  DIA  was  identified  in  the  VX  vapor,  as  shown  in  Figure  8,  at  a  nearly  constant 
rate  of  13.5  pg-min  ’,  with  an  integrated  mass  of  0.9  mg.  Assuming  the  results  from  the  GC-MS 
analysis  refleeted  the  molar  concentration  of  DIA,  and  that  an  insignificant  fraction  of  its 
starting  mass  was  consumed  during  the  first  day’s  experiment,  the  Rpr<;<i  was  42  pg-min  l.  The 
total  for  aeetone,  D1CDI,  and  DIA  together  was  4. 1  mg,  whieh  contrasted  with  the  ehange  in 
mass  of  6.9  mg  from  the  saturator  eell  weight.  The  Rpred  of  VX  was  0.20  pg-min  ',  and  the  fipre<j 
of  of  DEMP  was  0.24  pg-min  ’.  The  predicted  total  mass  of  both  compounds  was  0.03  mg. 

T  hroughout  the  series  of  experiments,  we  looked  for  evidence  of  the  2-(diiso- 
propylamino)ethanethiol  (VX  thiol)  in  the  VX  vapor  by  considering  the  following  eharaeteristies 
of  the  compound: 

•  Its  reported  area  peak  fraction  (from  GC-MS)  and  pure  vapor  pressure 
were  higher  than  those  of  DEMP  in  the  sample. 

•  At  15  °C,  and  with  an  assumed  mole  fraction  of  0.3 1%,  the  /?prca  of  VX 
thiol  was  0.30  pg-min  '  and  total  mass  consumed  would  have  been 
expeeted  to  have  been  0.03  mg. 

•  Within  the  fingerprint  region  of  VX,  the  most  intense  band  in  the 
absorptivity  eoeffieient  of  VX  thiol  was  observed  at  1 166  em  ',  with 

a  =  0.00042  (pmol/mol)  'm  ',  giving  a  prediction  of  A-  0.0002,  whieh 
was  just  equal  to  the  RMS  noise  level. 
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•  The  absorption  feature  of  VX  thiol  was  in  a  crowded  region  of  the 
spectrum  of  VX  vapor;  VX  has  a  stronger  band  at  1 160  cm 

•  D1CDI  has  a  band  at  1 169  em  1  that  is  nearly  as  intense  as  VX  thiol,  with 
a  weaker  feature  from  DEMP  also  nearby. 

•  A  slightly  weaker  absorption  feature  was  observed  for  VX  thiol  at 

1210  em  1  [a  =  0.00042  (pmol/mol)  'm  ']  that  may  lie  within  the  shoulder 
of  the  bands  associated  w  ith  the  P-0  stretch. 


Figure  8.  Spectrum  of  effluent  from  VX-filled  saturator  eell  acquired  on  26  February  2009.  The 
reetangle  highlights  an  absorption  band  near  700  cm  \  which  is  associated  w  ith  the  DIA  found  in 
the  vapor.  The  DICD1  peak  is  at  640  em  '. 

Of  the  remaining  experiments,  all  were  eondueted  with  the  saturator  cell  set  to 
23  °C.  Among  these,  the  next  eight  experiments  (encompassing  4—6,  10-12,  and  16-17  March 
2009)  were  run  with  the  earner  at  100  seem  and  the  diluent  at  2500  seem. 

3.5  Quantitation  of  Aeetone  by  Vapor-Phase  1R 

Among  the  ehemieals  found  in  VX  that  were  not  CW  agents,  the  behavior  of 
acetone  proved  to  be  the  most  interesting.  Initially,  the  mass  rate  of  aeetone  deeayed  rapidly  at 
23  °C  and  appeared  to  follow  the  expected  logarithmic  behavior.  Figure  9  shows  the 
experimental  data  from  4  to  5  Mareh  2009,  with  R  beginning  at  3 1 .4  pgmin  1  at  /  =  0  and 
declining  to  <1 .5  pg-min  1  at  t  =  125  min.  The  apparent  scattering  of  the  data  points  toward  the 
end  of  the  plot  arises  from  the  faet  that  R  was  calculated  from  signals  with  relatively  low  signal* 
to-noise  ratios  (SNRs).  The  final  data  point  shown  in  the  plot  was  calculated  from  an  \A  in 
which  /lmax  =  0.0008,  in  a  region  of  the  spectrum  with  an  RMS  noise  level  of  0.0001.  As  the 
eoneentration  of  aeetone  in  the  vapor  declined,  we  also  noted  the  presence  of  a  peak  centered 
around  1676  em  whieh  was  near  the  lower  frequency  shoulder  of  the  aeetone  earbonyl  stretch 
that  was  not  associated  with  an  aeetone  vibration.  The  frequency  of  the  vibration  would  be 
consistent  with  an  amide.  The  only  compound  in  this  class  that  was  reported  to  be  found  in  the 
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VX  was  1 ,3 -diisoprop ylurea.  Based  upon  GC-1R  spectra  of  the  compound  acquired  in  our 
laboratory,  1,3-diisopropylurea  likely  has  a  vapor-phase  amide  carbonyl  vibration  that  is  close  to 
1 700  cm  1  and  a  more  intense  band  near  1 503  cm  . 


Figure  9.  Mass  rate  of  acetone  observed  in  effluent  from  VX-filled  saturator  cell  on 
4  and  5  March  2009.  The  decay  in  the  output  of  the  acetone  appeared  to  follow  the  logarithmic 
trend  predicted  in  eq  5,  with  ln(7?)  =  at  +  b ,  where  a  =  -0.02525  and  b  =  3.446. 

On  6  March  2009,  the  output  of  acetone  observ  ed  in  effluent  from  the  VX-filled 
saturator  cell  exhibited  a  recovery  that  started  at  3.4  pgmin  Because  an  oxidation  reaction  to 
produce  the  acetone  appeared  to  be  unlikely,  the  most  likely  explanation  for  the  rebound  in  the 
compound  was  the  mass  transport  phenomena;  i.e.,  the  acetone  was  depleted  at  the  liquid-carrier 
interface  more  quickly  than  it  was  replaced  by  mixing  within  the  saturator  cell.  A  further 
recovery  of  acetone  was  expected  during  the  4  day  delay  in  the  experiments  from  6  to  10  March 
2009.  There  was  a  small  rebound  in  acetone  output,  from  just  under  0.4  pgmin  1  at  the  end  of 
6  March  to  0.5  pgmin  1  at  the  start  of  the  experiment  on  10  March.  The  decay  in  the  acetone 
continued  on  10  and  1 1  March,  and  by  12  March  the  compound  could  no  longer  be  detected  in 
the  vapor.  The  experiments  were  then  delayed  for  4  days.  Although  we  expected  that  a  small 
rebound  in  acetone  output  might  be  observed,  when  the  sample  was  reanalyzed  on  16  March,  we 
were  intrigued  to  note  a  dramatic  recovery  in  the  concentration  of  acetone  found  in  the  vapor 
from  the  saturator  cell.  This  is  illustrated  in  Figure  10,  which  shows  that  the  output  of  acetone 
was  nearly  as  high  on  16  March  as  it  was  on  4  March.  We  considered  the  possibility  that  the 
absorption  feature  could  have  come  from  another  similar  compound;  therefore,  we  verified  the 
identity  of  the  acetone  by  matching  other  absorption  features  to  the  reference  spectrum,  as  shown 
in  Figure  1 1.  Between  experiments,  the  saturator  cell  was  stored  in  a  laboratory  fume  hood,  with 
the  inlet  and  outlet  ports  sealed  with  0.25  in.  Swagelok  fittings  and  maintained  at  room 
temperature,  which  did  not  exceed  a  range  of  23  ±  3  °C. 
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At  the  conclusion  of  the  experiment  on  19  March,  low  concentrations  of  acetone 
were  still  detectable  in  the  VX  vapor  (albeit  <1  pg-min  ').  The  total  integrated  mass  of  acetone 
from  the  beginning  to  the  end  of  the  experiments  was  6.95  mg  or  0.000120  mol.  The  computed 
mass  fraction  was  0.092%  with  a  molar  fraction  of  0.45%.  At  a  mole  fraction  of  0.45%,  the 
Raoulf  s  law  prediction  yielded  an  initial  mass  rate  of  flo/Raouit  =  64.3  pg-min  1  at  5  °C  (the 
temperature  on  the  first  day  of  experimentation),  compared  with  an  experimental  value  of  Roitxp 
=  50.7  pg-min  1  (mass  rate  at  t  =  0  that  was  calculated  from  the  equation  in  Figure  7). 

Among  the  compounds  studied,  acetone  was  the  only  impurity  for  which 
independent  confirmation  of  its  molar  fraction  was  available  (through  NMR  spectroscopy). 
Although  GC-MS  offers  several  techniques  that  can  provide  more  rigorous  quantitation  of  trace 
compounds  in  a  solution  (through  the  use  of  external  and  internal  standards,  for  example),  such 
procedures  are  dependent  upon  the  compounds  and  matrix  and  may  require  extensive  methods 
development.  Such  an  extended  effort  involving  multiple  compounds  would  have  been  beyond 
the  scope  of  this  study.  Nevertheless,  the  availability  of  the  NMR  data,  at  least  for  acetone, 
proved  extremely  valuable  for  validation  of  the  results  from  1R  spectroscopy. 


Figure  10.  Mass  rate  of  acetone  in  effluent  from  VX-filled  saturator  cell  on  16  March  2009. 
Compared  with  the  12  March  results,  the  acetone  output  exhibited  a  dramatic  rebound. 


Figure 
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Spectrum  of  VX  effluent  and  acetone  reference  spectrum  showing  a  match. 


3.6  Quantitation  of  D1A  by  Vapor-Phase  1R 

As  noted  in  Table  1,  DIA  was  detected  by  GC-MS  analysis  in  the  VX  sample  at  a 
peak  area  fraction  of  0.3 1%.  If  it  is  assumed  that  this  value  more  closely  approximates  the  molar 
fraction  of  the  compound  in  the  VX,  the  predicted  initial  mass  rate  at  5  °C  becomes  R  — 

24.9  pg  min  as  noted  in  Section  3.4.  Although  eq  4  predicted  that  >1.5  mg  of  DIA  should 
have  evaporated  from  the  solution  during  the  first  day’s  experiment,  the  weakness  of  the 
absorption  feature  precluded  the  possibility  of  detecting  the  compound  at  that  temperature. 


DIA  was  first  seen  in  the  VX  vapor  at  1 5  °C.  The  SNR  was  only  5.  Nevertheless, 
the  availability  of  16  data  points,  the  statistical  advantage  of  signal  averaging  across  a  spectral 
region,  and  the  nearly  steady  rate  of  the  compound  enabled  us  to  calculate  a  mean  rate  of  13.5  ± 
2.4  gg  min  1  and  a  total  mass  of  0.87  mg.  This  was  significantly  less  than  the  predicted  value  of 
R  -  42.2  gg  min  1  (assuming  that  an  insignificant  mass  of  DIA  had  been  consumed  on  the  first 
day). 
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Figure  12.  Mass  rate  of  DIA  observed  in  effluent  from  VX-filled  saturator  eel!  from  4  to  6 
March  2009.  The  decay  in  the  output  of  the  compound  appeared  to  follow  the  logarithmic  trend 
predicted  in  eq  5,  with  ln(/?)  =  b0(  +  b\,  where  b0  =  -0.0063 1 1  and  b\  =  4.21 1. 

During  the  third  day’s  experiment,  with  the  temperature  of  the  saturator  cell  set  to 
23  °C,  the  output  of  DIA  began  to  follow  the  expected  logarithmic  decay,  as  shown  in  Figure  12 
which  contains  a  plot  of/?  for  the  three  experiments  run  between  4-6  March  2009.  The  decay 
equation  indicates  an  Ro  =  67.4  pg-rnin-1,  which  was  five  times  the  rate  obtained  at  15  °C.  The 
integrated  mass  of  DIA  was  8.8  mg.  Raoult’s  law  predicted  an  increase  proportional  to  the 
change  in  the  vapor  pressure  of  the  pure  compound,  9772  Pa/6613  Pa  =  1 .47,  which  was  more 
than  three  times  greater  than  the  A R  observed.  A  rebound  effect  similar  to  that  observed  w  ith  the 
acetone  was  possible,  and  the  grouped  data  points  from  the  experiments  run  on  the  3  days 
appear  to  reflect  a  tendency  for  the  rate  of  the  compound  to  “recover”  slightly  at  the  beginning  of 
the  operations  on  5  and  6  March.  However,  nitrogen-containing  compounds,  particularly  at  low 
concentrations  and  carrier  flows,  also  tend  to  be  “sticky”  in  our  system  and  required  significantly 
more  time  to  passivate  than  was  indicated  by  the  mixing  rate.  It  is  more  likely,  therefore,  that  the 
data  from  the  operations  at  5  and  1 5  °C  underestimated  the  mass  of  DIA.  Starting  on  6  March, 
the  total  time  of  each  run  was  increased,  which  gave  the  more  difficult  compounds  additional 
time  to  equilibrate  in  the  White  cell. 
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Figure  13.  Mass  rate  of  D1CDI  observed  in  effluent  from  VX-filled  saturator  cell  from  4  to  16 
March  2009.  The  decay  in  the  output  of  the  compound  appeared  to  follow  the  logarithmic  trend 
predicted  in  eq  5,  with  ln(ft)  =  b0t  +  b\,  where  /»o  —  —0.00 1018  and  b\  -  4.0358  and  t\n  -  680.7 
min. 


A  more  evident  rebound  effect  was  observed  when  the  mass  rate  of  D1A 
approximately  doubled  between  the  end  of  the  operation  on  6  March  and  the  start  of  the 
experiment  on  10  March  2009.  The  integrated  mass  of  the  D1A  from  the  two  experiments 
performed  on  10  and  1 1  March  was  3.8  mg.  After  1 1  March,  the  compound  was  no  longer 
detectable  in  the  VX  vapor. 

The  total  observed  mass  of  D1A  from  all  experiments  was  13.5  mg,  which  was 
equivalent  to  a  mass  fraction  of  0. 1 8%.  Based  upon  its  retention  time,  the  unknown  compound 
in  Tabic  1  (RT  =  20.37  min)  was  probably  bis[S-(2-diisopropyIamino)ethyl]  methyl 
phosphonodithiolate  (VX  bis).  Assuming  the  peak  area  fractions  in  Table  1  more  closely 
a  approximate  molar  fractions,  the  mean  molecular  weight  of  the  VX  sample  was  266.4.  The 

calculated  starting  molar  fraction  of  the  DLA  was  then  0.47%,  which  gave  a  predicted  value  of  R 
=  37.9  pgmin  1  at  the  beginning  of  the  experimental  series  (7*=  5  °C).  Although  a  minimum 
detection  limit  for  D1A  was  not  calculated,  we  believe  it  to  have  been  <13  pgmin  '.  The  failure 
to  detect  this  compound  in  the  VX  vapor  at  a  bath  temperature  of  5  °C  indicates  that  the  partial 
pressure  and  mass  rate  may  have  been  suppressed  in  the  solution  more  than  a  simple  Raoult’s 
law  calculation  would  have  predicted.  Even  at  23  °C,  the  initial  mass  rate  of  the  DIA  was  less 
than  two  times  that  predicted  at  a  temperature  of  5  °C. 
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3.7 


Quantitation  of  D1CD1  by  Vapor-Phase  I R 


When  acquiring  the  vapor-phase  absorptivity  coefficient  of  DICD1,  we  found  that, 
similar  to  the  testing  of  other  nitrogen  compounds,  D1CD1  required  more  time  than  other  less 
“sticky”  compounds  to  equilibrate  in  our  system.  It  did,  however,  have  a  strong  [a  =  0.00805 
(jimol/mol)  'm  '],  narrow  absorption  feature  near  2130  cm  1  in  an  uncrowded  region  of  the  1R 
spectrum.  Those  properties  made  D1CD1  a  good  compound  for  quantitation  using  1R 
spectroscopy,  and  allowed  us  to  compute  the  concentration  of  the  compound  for  the  full  series  of 
spectra  using  the  integrated  absorbance  across  the  band.  Because  of  the  relatively  high 
concentration  compared  with  the  other  volatile  compounds  (peak  area  fraction  by  GC-MS  = 
0.84%),  D1CD1  was  observed  at  a  high  SNR  ratio  (183),  even  at  the  end  of  the  series  of 
experiments.  In  fact,  evaluation  of  the  data  from  the  first  day’s  experiment  at  a  saturator  cell 
temperature  of  5  °C  indicated  that  the  starting  molar  concentration  of  the  compound  in  the  VX 
may  have  been  >2%. 

As  noted  in  Section  3.4,  the  mass  of  D1CD1  (computed  from  the  vapor-phase  IR 
spectra)  totaled  2.4  mg  during  the  first  two  experiments  at  5  and  15  °C. 

The  experiments  conducted  from  4  to  16  March  2009  were  all  run  with  a  carrier 
rate  of  100  seem  and  a  saturator  cell  temperature  of  23  °C.  In  accordance  with  eq  5,  we  expected 
to  see  a  logarithmic  decay  in  the  mass  rate  of  the  compound,  similar  to  the  behavior  observed 
with  the  acetone  and  DIA.  However,  contrary  to  our  expectation,  when  each  day’s  experiments 
were  plotted  separately,  they  appeared  initially  to  follow  a  linear  behavior.  Although  this  linear 
behavior  seemed  unlikely,  the  results  from  each  day’s  experiments  were  computed  individually, 
resulting  in  a  total  mass  of  35.8  mg  for  the  seven  experiments.  After  the  final  experiment  was 
completed  at  a  earner  rate  of  100  seem,  the  pooled  data  appeared  to  exhibit  the  predicted 
logarithmic  trend,  as  shown  in  Figure  13  (Section  3.6).  Using  the  equation  of  the  fitted  line  in 

J-997 

/(/)  =  35.4  mg,  little  had  changed  from  the 

0 

mass  that  was  obtained  by  computing  the  experiments  individually.  The  half-life  of  the 
compound  under  the  experimental  conditions  was  680.7  min.  The  hook-shaped  appearance  in 
the  early  part  of  the  data  from  each  day’s  run  reflects  the  additional  equilibration  time  required 
for  the  D1CD1,  which  has  two  nitrogen  atoms  per  molecule. 


For  the  final  set  of  three  experiments  beginning  on  17  March  2009,  carrier  flow 
rates  >100  seem  were  used.  As  described  in  Section  2.1,  higher  flow  rates  reduce  the  efficiency 
of  the  saturator  cell.  In  such  cases,  eq  4  must  be  rew  ritten  to  account  for  the  reduced  efficiency 
of  the  saturator  cell  at  the  higher  flow  rate 


n 


a 


PXn  . ,  _ 

a  a  earner  r 


(H) 


where  E  is  the  efficiency  of  the  saturator  cell  (from  Figure  2).  Because  eq  1 1  encompasses  a 
dynamic  environment,  it  is  used  to  predict  the  instantaneous  mass  rate  of  a  compound  only.  As 


30 


we  have  observed  through  working  with  other  compounds  at  high  carrier  rates,  the  decay  of  a 
volatile  impurity  in  a  compound  does  not  necessarily  follow  a  first-order  trend  with  respect  to  the 
natural  logarithm  of  its  mass  rate.  As  the  carrier  rate  increases,  the  phenomenology  of  mass 
transport  from  the  bulk  liquid  to  the  wick  and  then  to  the  vapor-phase  presumably  becomes  more 
complex,  leading  to  the  depletion  of  the  more  volatile  compounds  at  the  interface  between  the 
liquid  and  vapor.  This  trend  is  depicted  in  Figure  14,  which  shows  a  graph  of  the  mass  rate  of 
D1CD1  on  1 7  March  2009,  obtained  at  a  carrier  rate  of  500  seem.  The  mass  rate  of  the 
compound  on  1 8  and  19  March  2009  followed  a  similar  trend. 
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Figure  14.  Mass  rate  of  DICD1  observed  in  effluent  from  VX-filled  saturator  cell  on  17  March 
2009.  The  decay  in  the  output  of  the  compound  followed  the  trend  ln(fl)  =  hot1  +  b\t  +  b2,  where 
b0  =  -0.00001 299,  b\  =  -0.006986,  and  b2  =  4.2991 .  Total  integrated  mass  was  8.5  mg. 

Table  4  lists  the  integrated  mass  of  D1CD1  for  all  experiments  performed  in  this 
study.  As  noted  previously,  the  absorption  band  for  the  compound,  even  at  the  end  of  the  final 
experiment,  was  still  at  an  SNR  >10.  The  57.4  mg  shown  in  the  last  column  of  the  table  docs 
not.  therefore,  represent  the  original  starting  mass  of  DICD1  in  the  VX.  A  better  estimate  of  the 
total  D1CDI  starting  mass,  found  in  the  sample  of  VX  in  the  saturator  cell,  could  be  obtained  by 
combining  (1 )  the  masses  from  the  first  two  experiments  on  25  and  26  February  2009,  and  (2) 
the  mass  obtained  by  extending  and  integrating  the  equation  for  Figure  13  out  to  10  half-lives  (at 
which  point  the  mass  rate  would  have  decayed  to  0.5 ro  =  0.00098  of  the  starting  rate).  This 
resulted  in  a  total  starting  mass  of  0.6  +  1.8  +  55.8  =  58.2  mg,  which  is  a  molar  fraction  of  1.6  %. 
At  this  concentration,  the  Raoult’s  law  prediction  of  the  partial  pressure  at  5  °C  was  2.2  Pa,  and 
the  predicted  initial  mass  rate  of  D1CD1  at  50  seem  carrier  (the  conditions  of  the  first  experiment) 
was  R0  =  6.3  pgmin  l,  compared  with  the  observed  rate  of  8.8  pg-rnin  '. 
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Table  4.  Integrated  Mass  of  D1CDI  from  VX  Experiments,  25  February  to  19  Mareh  2009. 


Date 

25  Feb 

26  Feb 

4-16  Mar 

17  Mar 

18  Mar 

19  Mar 

Total 

Mass  (mg) 

0.6 

1.8 

35.4 

8.5 

7.1 

4.0 

57.4 

3.8  Quantitation  of  DEMP  and  VX  by  Vapor-Phase  1R 

As  noted  earlier,  neither  DEMP  nor  VX  was  observed  speetroseopieally  at 
saturator  eell  temperatures  of  5  or  15  °C.  The  conditions  used  in  the  experiments  run  from  4  to 
16  Mareh  2009  included  a  saturator  cell  temperature  of  23  °C  with  a  carrier  rate  of  100  seem  and 
a  diluent  rate  of  2500  seem.  The  predieted  mass  rates  were  1  and  0.3  pgmin  1  for  VX  and 
DEMP,  respectively.  Even  though  both  compounds  have  strong  absorption  features  near 
1040  em  ',  sueh  low  mass  rates  would  be  expected  to  result  in  a  value  of  A  =  0.0002  for  the  two 
compounds  combined.  Furthermore,  DIA  has  a  partially  overlapping  band  with  a  maximum  at 
1020  em  '.  At  an  RMS  noise  level  of  0.0001  w'ithin  the  speetral  region,  sueh  a  low  signal  would 
have  been  expeeted  to  be  difficult  to  deteet.  This  is  refleeted  in  Figure  15,  which  illustrates  the 
challenge  of  confirming  the  presence  of  the  two  phosphorous  compounds. 
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Figure  15.  Speetra  of  vapor  from  VX-filled  saturator  eell  (acquired  on  4  March  2009),  DEMP 
and  DIA.  At  a  predicted  value  of  A  =  0.0002  for  the  most  intense  absorption  band,  the  presence 
of  the  two  phosphorous  compounds  would  have  been  difficult  to  confirm  speetroseopieally. 

As  the  earner  rate  was  increased  to  500  seem  and  then  to  1000  seem,  the  presence 
of  DEMP  was  clearly  detectable  on  the  basis  of  the  P=0  streteh  at  1269  em  1  and  the  P-O-C 
streteh  in  the  1000  em  1  region.  Figure  16  shows  a  comparison  of  the  P-O-C  stretching  bands  in 
spectra  acquired  on  sueeessive  days  at  earrier  rates  of  100,  500,  and  1000  seem.  Integrated 
masses  for  the  experiments  run  with  earrier  rates  of  500  and  1000  seem  were  0.3  and  1.0  mg  on 
17  and  18  Mareh,  respectively.  Based  upon  the  appearanee  of  the  bands  shown  in  Figure  16  and 
the  lack  of  a  P~0  band  from  VX,  the  bands  shown  in  the  figure  likely  resulted  from  DEMP. 
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Figure  16.  Spectra  of  vapor  from  VX-filled  saturator  cell  at  carrier  rates  of  100,  500,  and 
1000  seem.  The  spectra,  which  were  acquired  at  30  min  into  each  experiment,  reflect  the 
increasing  concentration  of  DEMP  at  the  higher  carrier  rates.  The  shape  of  the  hands  in  this 
region,  as  well  as  the  absence  of  a  P=0  band  from  VX  (not  shown  above),  indicate  that  the 
features  likely  do  not  show  the  presence  of  VX  in  the  spectra. 

The  final  experiment  was  run  at  a  carrier  and  diluent  rate  of  2000  seem.  By  using 
the  high  carrier  rate,  we  hoped  to  increase  the  evaporation  of  VX  from  the  saturator  cell  and 
reduce  the  equilibration  time  in  our  system  by  providing  a  more  efficient  sweep  of  the  transfer 
lines  and  White  cell.  Analysis  of  the  spectra  taken  at  30  and  180  min  into  the  experiment 
(Figure  17)  indicate  that  the  high  flow  rates  effectively  increased  the  downstream  concentration 
of  VX  and  permitted  the  detection  and  identification  of  the  compound  in  the  spectra.  The 
identification  of  VX  was  facilitated  by  (1 )  the  presence  of  the  P=0  absorbance  bands  from  both 
DEMP  (1270  cm  ')  and  VX  (1250  cm  ')  at  frequencies  with  little  overlap,  and  (2)  the  change  in 
appearance  of  the  bands  in  the  1000  cm  1  region.  The  DEMP  band  has  a  shoulder  at  1070  cm 
in  contrast  to  VX,  which  has  only  the  single  band  at  1045  cm  *.  Even  at  the  higher  flow  rate,  VX 
was  still  a  “sticky”  compound,  and  the  increase  in  the  intensities  of  the  P-O-C  bands  between  30 
and  180  min  can  easily  be  seen,  along  with  a  further  change  in  their  overall  shape  to  a  more 
VX-like  appearance.  The  arrows  shown  in  Figure  17a  and  b  indicate  the  presence  of  an 
absorption  feature  in  the  reference  spectrum  near  1 228  cm  1  that  was  not  seen  in  the  spectra  of 
the  vapor  from  the  saturator  cell.  Evaluation  of  the  GC-1R  (vapor-phase)  spectra  of  VX, 
acquired  in  our  laboratory,  (Figure  18)  have  indicated  that  the  lower  frequency  band,  if  present  in 
the  vapor-phase  spectrum  of  the  compound,  is  much  weaker  than  that  seen  in  the  reference 
spectrum  (Figure  17).  The  lower  frequency  band  in  the  VX  reference  spectrum,  which  matches 
the  frequency  of  the  P=0  band  in  liquid-phase  VX,  likely  arose  from  aerosolized  VX  in  the 
vapor  used  to  generate  the  spectrum.  Furthermore,  the  observed  A  v  of  =20  cm  1  in  the 
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frequency  of  the  P=0  hand  in  the  vapor  compared  with  the  liquid-phase  spectrum  of  an 
alkylphosphonate  is  typical  of  other  compounds.' ' 
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Figure  17.  Spectra  of  vapor  from  weapon s-grade  VX  in  a  saturator  cell  at  (a)  30  min  and  (h) 

1 80  min.  The  figures  include  the  vapor-phase  reference  spectra  of  DEMP  (ECBC)  and  VX 
(PNNL).  The  arrows  highlight  an  absorption  feature  in  the  VX  reference  spectrum,  apparently 
associated  with  aerosolized  VX,  that  was  absent  in  the  spectrum  of  the  material  generated  using  a 
saturator  cell. 
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Figure  18.  IR  spectra  of  VX.  The  P=0  stretching  band  in  the  GC-IR  (vapor-phase)  spectrum 
(blue  band)  was  near  1250  cm  '.  This  contrasts  with  the  liquid-phase  spectrum  (black  band),  in 
which  the  band  was  near  1228  cm  '.  The  vapor-phase  absorptivity  coefficient  of  VX  from  the 
PNNL  database  (red  band)  appears  to  exhibit  P=0  absorption  features  from  both  liquid  and 
vapor-phase  VX. 
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The  predicted  values  for  the  mass  rate  and  absorbance  of  the  two  compounds  in 
the  final  experiment  arc  shown  in  Table  5.  Because  the  vapor  pressures  of  the  compounds  have 
been  well  established  by  measurements  conducted  at  ECBC,2,1  the  prediction  uncertainty  for  R 
is  primarily  dependent  upon  (1)  deviations  from  Raoult’s  law',  which  must  be  determined 
empirically,  and  (2)  uncertainties  in  the  concentrations  of  DEMP  and  VX.  We  can  assume  that 
the  purity  of  the  VX  is  probably  at  least  close  to  the  value  reported  by  using  the  GC-MS.  For 
DEMP,  however,  the  relationship  between  the  area  fraction  reported  by  using  the  GC-MS  and 
the  mass  or  mole  fraction  is  unknown.  Under  the  assumption  that  the  fragmentation  of  VX  and 
DEMP  by  the  electron  beam  and  the  chromatographic  behavior  of  the  two  compounds  were 
similar,  we  used  the  peak  area  fraction  from  the  GC-MS  analysis  to  compute  a  mass  fraction  for 
the  DEMP. 


The  vapor-phase  absorptivity  coefficient  of  DEMP  has  been  measured  by  the 
personnel  at  PNNL14  and  ECBC.1'  Within  the  region  that  includes  the  1045  cm  1  band, 
deviations  between  the  data  from  the  two  laboratories  was  only  2.7%.  Therefore,  the  uncertainty 
in  the  a  value  of  the  compound  has  little  effect  on  the  uncertainty  of  the  predicted  value  of  A. 
There  has  been  no  interlaboratory  comparison  of  the  vapor-phase  absorptivity  coefficient  of  VX. 
The  metadata  file  associated  with  the  data  contains  uncertainties  of:  Type  A  =  9.5%  and  Type  B 
<10%.  A  note  in  the  file  indicates  that,  “some  liquid  probably  condensed  on  White  cell  mirrors.” 
As  we  have  further  noted,  this  appears  to  be  confirmed  by  the  appearance  of  the  spectra  within 
the  fingerprint  region.  The  Type  B  uncertainty  that  was  reported  may,  therefore,  be  an  optimistic 
finding.  In  addition  to  the  DEMP  and  VX,  VX  thiol  has  a  maximum  absorption  band  at  1069 
cm  1  w  ith  a  shoulder  that  extends  beyond  1045  cm  1  and  likely  contributes  at  least  weakly  to  the 
spectra  in  Figure  17.  Given  the  associated  uncertainties,  the  prediction  values  in  the  last  column 
of  Table  5  are  remarkably  similar  to  the  actual  spectra. 


Table  5.  Predicted  Mass  Rates  and  Absorbance  Values  in  Vapor-Phase  Spectra  of  Weapons- 
Grade  VX.  The  values  shown  in  the  table  arc  for  the  experiment  conducted  at  a  carrier  rate  of 
2000  seem.  The  values  reflect  the  intensities  of  the  band  at  1045  cm  1  that  should  be  contributed 
by  DEMP  and  VX  at  the  predicted  mass  rates  shown. 


Compound 

R  (pg  min  ') 

A  (login) 

DEMP 

4.4 

0.0012 

VX 

15.0 

0.0018 

The  integrated  mass  of  DEMP  from  the  experiments  in  which  the  presence  of  the 
compound  could  be  identified  and  quantified  (17-19  March  2009)  was  2.6  mg.  On  the  basis  of 
the  relationship  between  the  P-O-C  and  P=0  bands,  and  the  contribution  to  the  P-O-C  band  from 
DEMP,  the  mass  of  VX  observed  on  19  March  2009  was  calculated  to  be  2.7  mg.  Assuming 
ideal  Raoult  behavior  and  a  starting  molar  fraction  of  0.1%  DEMP,  the  computed  mass  of  DEMP 
from  the  experiments  run  25  February  to  16  March  was  1.5  mg,  w'ith  a  total  of  4.0  mg  for  all 
experiments. 
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3.9 


Mass  Balance 


Table  6  provides  the  computed  mass  balance  for  all  compounds  having  available 
spectral  data.  The  total  shown  for  D1CD1  reflects  the  integrated  masses  from  the  first  two 
experiments  and  the  calculated  mass  of  the  compound,  assuming  the  equation  in  Figure  13  was 
used  to  integrate  to  10  half-lives.  The  masses  of  DEMP  and  VX  include  both  integrated  spectral 
data  and  assumptions  from  Raoult's  law  for  experiments  in  which  the  spectral  signatures  were 
too  weak  to  be  detected.  The  total  mass  of  the  compounds  shown  in  Table  6  was  88.2  mg,  \ 

compared  with  a  gravimetrically  determined  change  in  mass  of  108.8  mg  (calculated  from  the 
difference  between  the  starting  weight  of  the  filled  saturator  cell  on  25  February  and  the  ending 
weight  on  1 9  March).  The  compounds  listed  in  Table  6  accounted  for  8 1 .0%  of  the  gravimetric 
mass  change  over  the  experimental  time. 


Table  6.  Mass  Balance  of  Compounds  Observed  in  the  Vapor  from  the  VX  Saturator  Cell  and 
Gravimetric  Data  (Obtained  from  the  Starting  and  Ending  Weight  of  the  Saturator  Cell). 


Compound 

Acetone 

D1A 

D1CDI3 

DEMPb 

vxc 

Total 

Gravimetric 

Change 

Mass  (mg) 

7.0 

13.5 

58.2 

4.0 

5.5 

88.2 

108.8 

“Assumes  dala  from  4  to  16  March  2009  extended  to  10  half-lives 

includes  inlegraled  speelral  data  from  17  to  19  March  2009  and  assumes  ideal  Raoult’s  law  behavior  from  25 
February  lo  1 6  March  2009 

‘Includes  spectral  dala  for  19  March  2009  and  assumes  ideal  Raoult’s  law  behavior  from  25  February  lo  18  March 
2009 


We  were  unable  to  identify  a  unique  absorption  band  of  sufficient  intensity  under 
the  experimental  conditions  for  quantitation  of  VX  thiol.  Therefore,  because  we  had  no  means  to 
estimate  its  deviation  from  ideal  Raoult’s  law  behavior,  we  did  not  compute  the  evaporation  of 
the  compound  from  the  saturator  cell.  Other  compounds  (listed  in  Table  1)  may  have  contributed 
to  the  vapor  from  the  VX,  but  were  not  considered  in  this  report  or  shown  in  Table  6  because,  to 
our  knowledge,  neither  the  quantitative  vapor-phase  reference  spectra  nor  the  vapor  pressure  data 
within  the  range  of  our  tests  were  available.  The  compounds  that  most  likely  contributed  at  least 
part  of  the  mass  in  the  VX  vapor  included  2-(diisopropylamino)cthyl  vinyl  sulfide,  1,3- 
diisopropylurea,  and  diethyl  methylpyrophosphonate  (VX  pyro). 


4.  CONCLUSIONS 

Wcapons-grade  VX  contains  stabilizers  and  impurities  with  pure  component 
vapor  pressures  that  vary  from  that  of  VX  by  more  than  6  orders  of  magnitude.  The  wide  range 
of  volatilities  means  that  partial  pressures  for  even  trace  impurities  may  be  many  times  higher 
than  that  of  the  CW  agent.  When  the  agent  is  vaporized  through  evaporation,  a  complex  mixture 
of  materials  results  and,  therefore,  the  spectral  signatures  are  dominated  by  the  more  volatile 
compounds  rather  than  by  the  VX.  If  the  concentrations,  identities,  and  pure  component  vapor 
pressures  of  the  impurities  in  VX  are  known,  and  if  reference  spectra  for  these  compounds  are 
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available,  it  is  possible  to  generate  an  artificial  spectrum  of  the  vapor  that  more  closely 
approximates  the  aetual  vapor-phase  spectral  signature  of  the  agent. 

Using  relatively  simple  mathematical  procedures  based  upon  Raoult’s  and  Beer's 
laws,  we  were  able  to  use  vapor-phase  IR  spectroscopy  to  obtain  the  mass  balanee  of  the 
compounds  that  evaporated  from  weapons-grade  VX.  We  were  also  able  to  determine  the 
starting  concentrations  of  several  impurities  in  the  agent  that  were  difficult  to  measure  by  other 
analytical  techniques. 
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ACRONYMS  AND  ABBREVIATIONS 


CAS 

Chemical  Abstracts  Service 

CW 

chemical  warfare 

DEMP 

diethyl  methylphosphonate 

DIA 

diisopropylamine 

>  DICDI 

1 ,3-diisopropylcarbodiimide 

ECBC 

U.S.  Army  Edgewood  Chemical  Biological  Center 

FT1R 

Fourier  transform  infrared  spectroscopy 

GC-IR 

gas  chromatograph  infrared  analyzer 

GC-MS 

gas  chromatograph  mass  spectrometer 

HeNe 

helium-neon 

HgCdTe 

mercury-cadmium  telluride 

IR 

infrared 

NIST 

National  Institute  of  Standards  and  Technology 

NMR 

nuclear  magnetic  resonance 

PNNL 

Pacific  Northwest  National  Laboratory 

rms 

root  mean  square 

SNR 

signal-to-noise  ratio 

VX  bis 

if 

bis[S-(2-diisopropylamino)ethyl]  methyl  phosphonodithiolate 

VX  thiol 

2-(diisopropylamino)ethanethiol 

VX  pyro 

diethyl  methylpyrophosphonate 

43 


